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HEATTRANSFER IN A FISSIONING URANIUM PLASMA REACTOR CAVITY 

by Albert F. Kascak 
Lewis Research Center 

SUMMARY 

The cavity heat -transfer processes in a fissioning uranium plasma test reactor were 
analyzed. The thermal neutron flux was limited to about 10 3 neutrons/(cm )(sec) or 
less, which means that the heat-source strength was limited to 20 kW/g of fissionable 
material. The problem is not in cooling the cavity wall, but in insulating the uranium 
plasma to keep its temperature above its local condensation temper atureji The problem 
is especially severe at the outer boundaries of the plasma region. 

Two schemes were investigated. The first scheme was to use mirrored cavity walls 
to reflect the thermal radiation back into the uranium plasma. The second scheme was 
to seed the transpirational wall coolant so as to make it opaque to thermal radiation, thus 
insulating the hot plasma from the cold wall. The analytical approach taken was to ap- 
proximate thermal radiation by a diffusion process with appropriate "jumped” boundary 
conditions on temperature at the edge of the plasma region 

For a mirrored metallic wall, a realistic wall reflectivity is about 90 percent. For 
this case the operating pressure was taken to be a few hundreds of atmospheres. At a 
neutron flux of 10^ neutrons/(cm^)(sec), the temperature at the edge of the plasma was 
calculated to be within 500 K of the condensation temperature. The corresponding cavity 
power was in the tens of megawatts. This operating region is, at best, marginal. If the 
mirrored cavity wall could be made with a reflectivity of over 95 percent, the required 
power could be reduced by at least a factor of 2 and the concept would become much more 
attractive. 

In the second scheme when seeded hydrogen was used as the wall coolant gas the 
edge temperature of the plasma was very close to the condensation temperature. When 
seeded argon was used as the wall coolant gas, the edge temperature of the plasma was 
much higher than condensation temperature. This result was explained with a simple 
one -dimensional model which indicated that increasing either the gas molecular weight or • 
the seed mass fraction increased the edge temperature of the uranium plasma. Heavily 
seeded argon used as the transpirational wall coolant appears to be the better scheme to 
keep the edge of the uranium plasma from condensing for a given reactor thermal neutron 
flux, pressure level, and test region size. 



INTRODUCTION 


The Fissioning Uranium Plasma Reactor Facility was a proposal to use a solid -core 
nuclear reactor to generate neutrons in a central cavity of the reactor (ref. 1). These 
neutrons produce fissions in a uranium plasma located in this cavity, shown in figure 1. 
The heat transfer from this plasma to a cavity coolant is the subject of this report. 

Like all reactors the power density for this reactor is equal to the neutron flux times 
the density of the fissionable material. Since the density of the uranium plasma is much 
less than the density of the solid core, the power generated in the uranium plasma is 
much lower than the power generated in the solid core. The power level of the solid core 
is limited; thus, the power level of the uranium plasma is also limited - but at a much 
lower value. The problem is not in cooling the cavity wall, but in insulating the uranium 
plasma to keep its temperature above its local condensation temperature. 

This can be done ir. oeverarways. First, the inside cavity wall can have a mirrored 
surface, which would reflect the radiation from the uranium plasma back into itself. 
Second, transpirationally cooling the walls with a gas seeded with small solid particles 
would insulate the hot plasma from the cold wall. (Solid particles are used to increase 
. the gas opacity to thermal radiation.) The objective of this study is to determine the re- 
quired wall reflectivity and the choice of wall coolant, the amount of solid seed, and/or 
the radial flow profile desired. 

The approach is to approximate thermal radiation by a diffusion process with appro- 
priate jump boundary conditions. The properties used for hydrogen and uranium are the 
same as those used in reference 2. The properties of argon are extrapolated from ref- 
erence 3. 


ANALYSIS 

A heat -transfer model of the cavity of the Fission Uranium Plasma Reactor Facility 
is shown in figure 2. Flow (either hydrogen or argon) enters a slotted wall from a 
higher pressure plenum at a temperature T . It then picks up the amount of heat de- 
posited on the wall, leaving the wall at a temperature T , and flowing into a lower 
pressure cavity. In terms of enthalpy h this heat balance is 

q w = -(pv) w (h w - h p ) (1) 


(The wall coolant is nonseeded in the mirrored wall case or seeded with depleted ura- 
nium in the other case. ) 
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Inside the cavity, there is an energy balance between the radially inward convection 
toward the plasma and the radially outward heat flux toward the wall. In vector notation 
this can be written as 

0 = (pv) . Vh + V • q - Q’” (2) 

where Q ,M is the internal heat-generation rate (or nuclear heat source). This heat 
source is proportional to neutron flux times the density of the fissionable material. If 
the neutron flux is not attenuated by the plasma, the internal heat generation from the 
nuclear fission is 


Q'” = hPf (3) 

where 77 is a constant with units of w/kg; and p.j is the density of fissionable material. 

The flow pv and the density of fissionable material are functions of radial position. 
Inside the plasma the flow is zero, the properties are just a function of temperature, and 
the heat flux can be written as a diffusion process. 


'-"(Tr-*-)" 


A variable g can be defined as 


■f( ! 


g- / I2£El + k dT 
3a / 


and then the heat flux becomes 


q = -Vg 


The energy equation (eq. (2)) becomes 


V‘g = -vp { 


In Cartesian coordinates this equation can be integrated and becomes 



= +27j(f^ - f) = (-q ) 2 


(4) 


(5) 


( 6 ) 


(7) 


( 8 a) 
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where f is defined by 


= / S 

Jo 


f ' /„ Pf dg 


(8b) 


Taking the negative square root of equation (8a) and substituting into equation (6) 
yields 


q = V277(f t - f) 


(9a) 


Integrating equation (9a) yields 



(9b) 


Evaluating equations (9) at the edge of the plasma and eliminating 7 ? gives the following 
relation for the centerline conditions: 


q e x e 



dg 

Vf t - f 


The power per unit mass of fissionable material is then 



2(f t - V 


(10a) 


(10b) 


If it is assumed that the nuclear heat source is a constant, the power per unit mass 
of fissionable material is 


where is defined by 



(11a) 
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(lib) 


P f 


1 



p f dx 


and the profile used to evaluate equation (lib) is 


g = g e 




(11c) 


A comparison of equations (10) and (11) gives an estimate of the error involved with 
the assumption of a constant nuclear heat source. 

If we assume a constant nuclear heat source and spherical geometry, the heat flux 
becomes 


q = q e 


r?p f R 


VV 3 


(12a) 


where p^ is defined by 


- 3 P e 


p f R 2 dR 


(12b) 


and the profile used to evaluate equation (12b) is 


= g - (R 2 - R 2 ) 

e 2R_ v e/ 


(12c) 


The edge condition can be calculated by using the technique of reference 4. 

,4 


4aaTg AR = 2a(q+ + q'^AR + Q*” AR 


(13a) 


and 


^e ^e " ^e 


(13b) 


Equations (13) can be combined and result in 
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(14) 


ctT 


4 

e 



If the wall coolant is transparent and the walls are specularly reflecting, the tech- 
nique of reference 5 gives 


q e 




V 




+ (1 " e w ^ 


+ 

w 


(15a) 


(15b) 


(15c) 


and 


q 


w 




(15d) 


where the radiation is partitioned between the inner uranium plasma and outer spherical 
wall, and between the outer spherical wall and itself. Equations (15) can be substituted 
into equation (14) and result in 


aT 


4 

e 



+ 




Q Mt 

4a 


(16) 


Equation (16) is also true for slab geometry. 

If the wall coolant is not transparent but very opaque, the heat flux is a diffusion 
process. The edge condition can then be approximated by assuming 

q e w ( 17a) 


and 


ctT 


4 

e 


CTT 0 + 2 


Q, 


e 

4a 


(17b) 
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where Tq is the coolant -side temperature of the interface. Equation (17a) is true if the 
gradient of the temperature, on the wall coolant side of the interface, is small. 

The nuclear heat source in the region outside the fissioning plasma is zero. Thus, 
for spherical geometry; equations (2) and (6) can be written as 

0 = pv^--i± (R 2 q) (18a) 

dR r 2 dR 


and 



dr 


where the properties are assumed to be functions of temperature but not position. Given 
the inside and outside wall temperature and with the use of equation (1), equations (18) 
can be solved as an initial -value problem from the wall inward. A jump boundary con- 
dition at the wall is not necessary since convection is dominant. 

To gain insight into this problem, assume the region between the plasma and the 
wall can be subdivided into two regions (see fig. 2(b)). Near the plasma, conduction is 
dominant and convection can be neglected. Near the wall, convection is dominant. As- 
sume the flow enters the cavity, flows inward at a constant rate through the first region, 
then kirns at the interface of the two regions, and flows out the exhaust nozzle of the 
cavity. 

Assume the transpirational wall coolant is a mixture of solid particles and a perfect 

transparent gas. The radius of the particles is R_ and the density is p_. Assume the 

s s 

gas is monatomic with a molecular weight equal to M and at a pressure of P and a 
temperature of T. The mass fraction y is the mass of solid particles per unit mass of 
coolant. The specific heat is 


C =if£Vl-y) (19a) 

p 2 \M/ 


where the specific heat of the solid particles has been neglected and 3? is the gas con- 
stant. The absorption coefficient is 


a = 




(19b) 


If molecular conduction is neglected, the conductivity becomes the radiative conductivity 
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(19c) 


k = 


64R s p s #crT 
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For a slab geometry, the heat flux in the conduction region is a constant. 



This can be integrated to give 


k T 
w w 



(20a) 


(20b) 


where the jump boundary condition is neglected and the subscript t refers to the fluid 
turning point at the interface. If there is no transpirational coolant, the turning point is 
the wall. 


q e (x w - x e> 

Vw 



(20c) 


The maximum heat flux at the wall, considering only transpirational cooling, is 

q w = -P vC p T w (21a) 

where the plenum temperature is assumed to be zero. When equation (21a) is used as a 
boundary condition, the heat flux in the convection region can be obtained by integrating 
equation (2). 

q = -pvC p T (21b) 

Evaluating equation (21b) at the turning point yields 

q e = -pvC pT t (21c) 

Using the conduction formulation of heat flux and equation (21b) yields 
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r 
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(24c) 


N c = -?(l - 4V)- 1 / 4 


where 


45pvPy(x - x ) 

N c = — 

128 Vs" T t 


(24d) 


For the case of no transpirational cooling, equation (24a) becomes 

i c = -(e 5 - 1) 

5 

For the case of a large amount of transpirational cooling, equations (24) become 


(25) 


« - «c + T (- N cf /4 


(26a) 


r = 




(26b) 


N c = 


-i 


V2 (- N^) 1 / 4 

Solving equations (26a) and (26b) for i; yields 


(26c) 


i = 


5 - r 


(27) 


The calculational procedure is summarized as follows: First, the assumption of con- 
stant heat generation in slab geometry was checked. The concept of a mirrored wall in 
spherical geometry was then investigated. The pros and cons of seeded argon versus 
seeded hydrogen were studied and explained with a simple one -dimensional model. 


DISCUSSION 

The Fissioning Uranium Plasma Reactor Facility was a proposal to use a solid-core 
nuclear reactor to ’’drive” a uranium plasma located in a central cavity (ref. 1 and 
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fig. 1). The operation of this fissioning uranium plasma was energy limited. The prob- 
lem addressed is not in cooling the cavity wall, but in insulating the uranium plasma to 
keep its temperature above its local condensation temperature. A heat -transfer model, 
shown in figure 2, was analyzed by using a diffusion approximation of thermal radiation 
with appropriate jump boundary conditions. The properties used were similar to those 
used in references 2 and 3. 

The calculational procedure was first to check the assumption of constant heat gen- 
eration in slab geometry. (Heat generation should be proportional to the density or to 
one over the temperature. ) This was done by comparing equations (10) (variable heat 
source) to equations (11) (constant heat source). The result of this calculation is shown 
in figure 3 with the solid lines being the variable-heat-source solution and the dashed 
lines being the constant-heat-source approximation. Over the pressures and the edge 
heat fluxes of interest, the heat-source- strength rj can be approximated by the constant- 
heat-source method. The same thing was noticed for a few spot calculations in spherical 
geometry; and therefore, the rest of the calculations use the constant- heat-source approx 
imation. 

To investigate the concept of a mirrored cavity wall, equations (12) and (16) were 
used to generate figure 4. Metallic walls are generally about 90 percent reflective. The 
wall temperature was set at 500 K. (Wall temperatures as high as 2000 K resulted in 
only slight changes from the 500 K nominal wall case. ) The uranium plasma was located 
in the central 25 percent of the cavity volume. The region between the plasma and the 
cavity wall was transparent. Conduction and convection off the plasma were neglected. 

Figure 4(a) shows the necessary heat-source strength to keep the edge of the ura- 
nium plasma 2000, 1000, 500, or 0 K above the condensation temperature. Reference 1 
indicates that about 10 0 neutrons/ (cm )(sec) is the upper limit on the thermal neutron 

17 

flux. This corresponds to about 2x10 W/kg, (20 kW/g), as the upper limit of the heat- 

n 

source strength. The only region in figure 4(a) which is less than 2x10 WAg is in the 
hundreds of atmospheres pressure, and the edge temperature of the uranium plasma 
must be within 500 K of the condensation temperature. From figures 4(b) and (c), for 
this region, the cavity power is between 5 to 20 MW and the wall heat flux is less than 

9 

2 kW/cm . This operational region is, at best, marginal. 

The performance of the mirrored cavity wall concept is a strong function of the wall 
reflectivity (or one minus the wall emissivity). Figure 5 shows the effect on the neces- 
sary heat-source strength of varying the wall emissivity. The calculation was done at 
100 atmospheres and with the edge temperature of the plasma at the condensation tem- 
perature. If the mirrored cavity wall could be made with a reflectivity of over 95 per- 
cent, the concept would become much more attractive. 

The alternative to using mirrored cavity walls is to use an opaque gas flowing ra- 
dially inward to convect the heat back toward the uranium plasma. Solid seed particles 
are used to increase the gas opacity to thermal radiation. There is a minimum gas flow 
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necessary for hydrodynamic stability of the uranium plasma or for wall cooling. The 
choice of wall coolant gas, the amount of solid seed, and the desired radial flow profile 
are investigated herein. 

Figure 6 shows some typical cavity temperature profiles for hydrogen and argon 
wall cooling gases. The gases were seeded to the same mass fraction, 0. 25, and had 
the same nondimensional flow profile (i. e. , a linear decrease to zero at the edge of the 
plasma). The flow through the wall was adjusted so that the hot gas stood off the wall 
about the same distance. For the hydrogen case, the edge temperature of the plasma 
was very close to the condensation temperature. For the argon case, the edge temper- 
ature of the plasma was much higher than the condensation temperature. 

Figure 7 shows the maximum uranium plasma edge temperature as a function of 
thermal neutron flux for hydrogen and argon for pressures of 133 and 200 atmospheres. 
(Thermal neutron flux is proportional to heat -source strength, where the constant of pro- 
portionality is higher for argon, since hydrogen absorbs some neutrons.) For hydrogen 
as the wall coolant, the edge temperature of the uranium plasma is very close to the 
condensation temperature. Argon is obviously the best choice of wall coolant from 
this viewpoint. 

To explain this result, the simple one -dimensional model (eqs. (24)) was used to 
generate figure 8. To minimize the reactor power for a given plasma edge temperature, 
wall temperature, seed fraction, wall coolant, and geometry, the product, -N c r, must 
be a minimum. Therefore, it is desired to use the least possible amount of flow through 
the wall and to have this flow turn as soon as possible. Obviously, N c and r are not in- 
dependent variables, as shown in figure 8. The minimum value of £ is £ . If the con- 
densation temperature is about 5000 K and the maximum wall temperature is about 1000 K, 
the value of e is greater than 5. Using equation (25) gives the value of ij c as greater 
than about 625. From figure 8, N c r must be greater than 100 before £ differs appre- 
ciably from £ . The maximum value of 4 is given from equation (27) when r is equal 
to 1 and becomes £ = 1. 25 £ c . Therefore, can only vary by 25 percent. 

Since £ is relatively constant, the minimum reactor power results from a minimum 
k w in equation (22). From equation (19c), k w can be reduced by increasing either the 
seed mass fraction or the gas molecular weight. Changing the molecular weight from 
that of hydrogen to argon results in a drop in reactor power of a factor of 15. Doubling 
the seed mass fraction from 0. 25 to 0. 50 results in a drop in reactor power of a factor 
of 3. 


CONCLUSIONS 

The cavity heat transfer in a fissioning uranium plasma reactor is analyzed in this 
report. Two schemes were investigated to keep the plasma edge temperature above 
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condensation. First was the mirrored -cavity -wall scheme, which wovdd reflect the ther- 
mal radiation back into the uranium plasma. Second was the seeded transpirational wall 
coolant, which would convect the heat back towards the plasma. 

The analysis resulted in the following conclusions for the mirrored-cavity-wall 
scheme: 

1. For a realistic wall reflectivity of 90 percent, the operating pressure was in the 

hundreds of atmospheres, the edge temperature of the plasma was within 500 K of the 

condensation temperature, the cavity power was in the tens of megawatts, and the wall 

2 

heat flux was less than 2 kW/cm . 

2. If the mirrored cavity wall could be made with a reflectivity of over 95 percent, 
the required power could be reduced by at least a factor of 2 and the concept would be- 
come much more attractive. 

The analysis resulted in the following conclusions for the seeded -transpirational - 
wall -coolant scheme: 

1. For seeded hydrogen as the wall coolant gas, the edge temperature of the plasma 
was very close to the condensation temperature. 

2. For seeded argon as the wall coolant gas, the edge temperature of the uranium 
plasma was much higher than the condensation temperature. 

3. Increasing either the gas molecular weight or the seed mass fraction increased 
the edge temperature of the uranium plasma. 

Seeded argon used as the transpirational wall coolant was the best scheme to keep 
the edge of the uranium plasma from condensing for a given reactor flux. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 26, 1973, 

503-04. 
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APPENDIX - SYMBOLS 


A area, 

a absorption coefficient, m~* 

C specific heat, J/(kg)(K) 

r 

f auxiliary function 

g auxiliary function 

h enthalpy , J /kg 

k thermal conductivity, J/(m)(sec)(K) 

m molecular weight 

N c nondimensional flow parameter, a measure of convection to conduction 

o 

P pressure, N/m 

3 

Q'" internal heat -generation rate, W/m 

q heat flux, W/m^ 

R radius, m 

St gas constant, 8. 3143 J/(K)(mol) 

r nondimensional depth of penetration of wall coolant flow 

T temperature, K 

v velocity, m/sec 

x distance, m 

y mass fraction of seed in wall coolant 

e ratio of plasma edge temperature to wall temperature 
emissivity of cavity wall 

r\ heat-source strength per unit mass of fissionable material, W/kg 

£ nondimensional heat flux 

£ nondimensional heat flux with no flow 

s c 

3 

p density, kg/m 

O 

p average density, kg/m 

a Stefan -Boltzmann constant, 5.6697x10"®, w/(m^)(K^) 
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Subscripts: 

t center of plasma 

e edge of plasma 

f fissionable material 

mol molecular 

p plenum 

s seed 

t turning point of wall coolant 
w cavity wall 

0 wall coolant adjacent to plasma 

Superscript: 

± component in plus or minus direction 
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Figure 3. - Comparison of variable heat source (proportional 
to density of fissionable material) to constant heat source 
as function of plasma edge heat flux for various pressures, 
for the reactor shown in figure 1. 


<b) Slab geometry. 

Figure Z - Heat-transfer model. 
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(a) Heat-source strength as function of pressure. 


(b) Cavity power as function of pressure. 
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<c) Cavity wall heat flux as function of pressure. 

gure 4. - Necessary heat-source strength, cavity power, or cavity wall heat flux for uranium plasma inside a 0. 61-meter (2-ft) diameter cavity with 90-percent 
reflecting walls and a transparent wall coolant gas to keep plasma edge temperature 2000. 1000, 500, or 0 K above the condensation temperature. 
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